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Chiral edge states are a hallmark of quantum Hall physics. In electronic systems, they appear as a 
macroscopic consequence of the cyclotron orbits induced by a magnetic field, which are naturally 
truncated at the physical boundary of the sample. Here we report on the experimental realization of 
chiral edge states in a ribbon geometry with an ultracold gas of neutral fermions subjected to an 
artificial gauge field. By imaging individual sites along a synthetic dimension, we detect the existence 
of the edge states, investigate the onset of chirality as a function of the bulk-edge coupling, and 
observe the edge-cyclotron orbits induced during a quench dynamics. The realization of fermionic 
chiral edge states is a fundamental achievement, which opens the door towards experiments including 
edge state interferometry and the study of non-Abelian anyons in atomic systems. 


Ultracold atoms in optical lattices represent an ideal 
platform to investigate the physics of condensed-matter 
problems in a fully tunable, controllable environment mm- 
One of the remarkable achievements in recent years has 
been the realization of synthetic background gauge fields, 
akin to magnetic fields in electronic systems. Indeed, by 
exploiting light-matter interaction, it is possible to imprint 
a Peierls phase onto the atomic wavefunction, which is 
analogous to the Aharanov-Bohm phase experienced by 
a charged particle in a magnetic field ISHg. These gauge 
fields, first synthesized in Bose-Einstein condensates [6], 
have recently allowed for the realization of the Harper- 
Hofstadter Hamiltonian in ultracold bosonic 2D lattice 
gases d ig, paving the way towards the investigation 
of different forms of bulk topological matter in bosonic 
atomic systems [319]. In the present work we are instead 
interested in the edge properties of fermionic systems 
under the effects of a synthetic gauge field. Fermionic edge 
states are a fundamental feature of 2D topological states of 
matter, such as quantum Hall and chiral spin liquids [ini 
m- Moreover, they are robust against changing the 
geometry of the system by keeping its topology, and can be 
observed even on Hall ribbons [12]. In addition, they offer 
very attractive perspectives in quantum science, such as 
the realization of robust quantum information buses m, 
and they are ideal starting points for the realization of 
non-Abelian anyons akin to Majorana fermions nans]. 

Here, we report the observation of chiral edge states 
in a system of neutral fermions subjected to a synthetic 
magnetic field. We exploit the high level of control in 
our system to investigate the emergence of chirality as a 
function of the Hamiltonian couplings. These results have 
been enabled by an innovative experimental approach, 
where an internal (nuclear spin) degree of freedom of the 
atoms is used to encode a lattice structure lying in an 
“extra dimension” [12], providing direct access to edge 
physics. In addition, we validate the chiral nature of our 



FIG. 1. A synthetic gauge field in a synthetic dimen¬ 
sion. a. We confine the motion of fermionic ultracold atoms 
in a hybrid lattice, generated by an optical lattice along a real 
direction x with tunneling t, and by a laser-induced hopping 
between spin states along a synthetic direction m. By inducing 
a complex tunneling along m, the atom wavefunction 

acquires a phase (f per plaquette, mimicking the effect of a 
transverse magnetic field B on effectively charged particles, b. 
Scheme of the ^^^Yb nuclear spin states and Raman transitions 
used in the experiment. 


system by performing quench dynamics, demonstrating 
how the particle motion shows edge-cyclotron orbits [T6] . 

We synthesize a system of fermionic particles in an 
atomic Hall ribbon of tunable width pierced by an effec¬ 
tive gauge field. One dimension is realized by an optical 
lattice, which induces a real tunneling t between differ¬ 
ent sites along direction x (see Fig. la). The different 
internal spin states are coupled by a two-photon Raman 
transition, which provides a coherent controllable cou¬ 
pling between different spin components. This can 

be interpreted as a complex tunneling amplitude between 
adjacent sites of an ”extra-dimensional” lattice Haim. 
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FIG. 2. Chiral dynamics in 2-leg ladders, a. The upper panel shows false-color time-of-flight images of the atoms in 
the m = —5/2 and m = —1/2 legs (averages of ~ 30 realizations). The central panels show the integrated lattice momentum 
distribution n{k) and the lower panels show h{k) (the numbers in the insets are the values of J determined from h{k)). 
Experimental parameters: VLi — 2 ti x 489 Hz, t — 271 x 134 Hz, VLi/t — 3.65, (p = 0.377r. b. Time-of-flight images and h{k) of 
the m = —1/2 leg for opposite values of the effective magnetic field. Experimental parameters: = 27r x 394 Hz, t = 27r x 87 

Hz, Qi/t — 4.53, ip — ±0.377r. c. Sketch of the 2-leg ladder configuration realized for this experiment. The arrows are a pictorial 
representation of the chiral currents, d. The circles show experimental values of \J\ for the m = —1/2 leg as a function of Qi/t 
(averages of datasets taken for ip = ±0.377r). The error bars are obtained with a bootstrapping method applied on ~ 30 different 
measurements. The shaded area is the result of a numerical simulation including thermal fluctuations (see Supplementary 
Materials). 


Eurthermore, the phase imprinting laid out by the Raman 
beams amounts to the synthesis of an effective magnetic 
field for effectively charged particles [4] with flux ip/27r 
(in units of the magnetic flux quantum) per plaquette (see 
Supplementary Materials). The Hamiltonian describing 
the system is 


H = 


EE[ 

j a 


EE 


2 


+ h.C.) + 


( 1 ) 


where c^j are fermionic creation (annihilation) op¬ 

erators on the site (j, a) in the real (j) and synthetic 
{a = 1,2,3) dimension, and n = The first line 

describes the dynamics along x, where t can be tuned 
by changing the intensity of the optical lattice beams. 
The dynamics along ifi is encoded in the second line: 
can be controlled by changing the power of the Raman 
beams, while ip can be tuned by changing their angle (see 
Supplementary Materials). Besides the tunneling terms, 
jij describes a weak trapping potential along x, while 
accounts for a state-dependent light shift, providing an 
energy offset along m. In our experiment we realize large 
synthetic magnetic fields corresponding to if 0.377r per 
plaquette. The carriers are constituted by alkaline-earth¬ 
like ^^^Yb atoms, initially prepared in a degenerate Eermi 


gas, and the sites of the synthetic dimension (see Eig. lb) 
are encoded in a subset of spin states {m} out of the 
1 = 5/2 nuclear spin manifold, thus providing up to 6-leg 
ladders. These atoms show SU(A/')-invariant interactions 
m, inhibiting the redistribution of the atoms among the 
different synthetic sites by collisional processes [iai2oi. 

The key advantage offered by the implementation of 
the lattice in a real+synthetic space is the possibility to 
work with a finite-sized system with sharp and addressable 
edges. Specifically, we focus on elementary configurations 
made up of fermionic ladders with a small number of 
“legs” connected by a tunnel coupling between them. A 
leg is constituted by a ID chain of atoms trapped in the 
sites of the real lattice in a specific spin state, whereas 
the “rungs” are provided by the synthetic tunneling (see 
Fig. la). The number of legs can be set by controlling 
the light shifts in such a way as to choose the number 
of spin states that are coupled by the Raman lasers (see 
Supplementary Materials). 

We first consider the case of a 2-leg ladder constituted 
by the nuclear spin states m = —5/2 and m = —1/2. A 
quantum degenerate ^^^Yb Fermi gas, at an initial tem¬ 
perature T cx 0.2 Tf (where is the Fermi temperature), 
is first spin-polarized in m = — 5/2. By slowly turning on 
the intensity of the real-space lattice, we prepare a system 
of ladders where all the atoms occupy the m = —5/2 leg 
with less than one atom per site (i.e. in a conductive 
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FIG. 3. Chiral edge currents in a 3-leg ladder, a. Experimental time-of-flight images (up), n{k) (center) and h(k) for 
each of the three legs m — —5/2, m = —1/2 and m = +3/2 constituting the ladder (the numbers in the insets of the graphs are 
the values of J determined from h{k)). Experimental parameters: = 27r x 620 Hz, t = 27r x 94 Hz, fli/t = 6.60, (p = 0.377r. 

b. Sketch of the 3-leg ladder configuration realized for this experiment, c. The circles show experimental values of J for each 
leg as a function of VLi/t. The shaded areas are the results of a numerical simulation (see Supplementary Materials). Eor 
both experimental and simulation data, the blue, green, and red color correspond to m = —5/2, m = —1/2, and m = +3/2 
respectively. 


metallic state). Then, by controlling the intensity and 
frequency of the Raman beams (see Supplementary Ma¬ 
terials), we slowly activate the tunnel coupling between 
the legs in such a way as to load the Hall ribbons in an 
equilibrium state, without populating excited bands. 

Despite the absence of a real bulk region, this 2-leg 
configuration is expected to support chiral currents with 
atoms flowing in opposite directions along the legs (see 
Fig. 2c), as investigated recently in bosonic systems [21]. 
In order to witness this property, we measure the rela¬ 
tive motion of the atoms in the two legs by spin-selective 
imaging of the lattice momentum distribution, obtained 
switching off the synthetic coupling and releasing the 
atoms from the lattice. In Fig. 2a (upper panel) we show 
two time-of-fiight images referring to both the m = —5/2 
and the m = —1/2 legs (see Fig. 2c) for = 27r x 489 
Hz and t = 27 t x 134 Hz (Oi/t = 3.65). Here we are 
interested only in direction +, which reflects the distri¬ 
bution of the lattice momenta k along the legs (in units 
of the real-lattice wavenumber ki, = ir/d, where d is the 
real-lattice spacing). The lattice momentum distribution 
in the transverse direction is a uniform square due to 
the presence of an additional optical lattice, required to 
suppress the dynamics along the y and z real directions 
(see Supplementary Materials). The central panel of Fig. 
2a shows the lattice momentum distribution n{k) after 
integration of the images along the transverse direction 


and normalization according to J n{k)dk = 1. We observe 
a clear asymmetry in n{k) (similarly to what reported 
in experiments with spin-orbit coupling in harmonically 
trapped gases [22H24]), which we characterize by defining 
the function 

h{k) = n{k) — n{—k) , (2) 

reported in the lower panel of Fig. 2a. The quantity 
J = /q h{k)dk provides a measurement of the lattice 
momentum unbalance and quantifies the strength of the 
chiral motion of the particles along the two legs. The 
values J = +0.056(3) for m = —5/2 and J = —0.060(7) 
for m = —1/2 are approximately equal in intensity and 
opposite in sign, directly evidencing the presence of chiral¬ 
ity in the system. We also perform the same experiment 
with a reversed direction of the synthetic magnetic field 
B (see Fig. 2b), observing a change of sign in J, corre¬ 
sponding to currents circulating in the opposite direction. 
This behavior confirms the interpretation of our data in 
terms of chiral currents induced by a synthetic magnetic 
field in a synthetic 2D lattice. 

The stability of chiral edge states in fermionic systems 
is of key importance, e.g. for quantum information ap¬ 
plications. In our system, the appearance of a chiral 
behavior is governed by several key parameters, including 
the ratio fli/t, the Fermi energy Ep and the flux p. Their 
ample tunability can be used to investigate the rise and 
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FIG. 4. Edge-cyclotron orbits, a. Time dependence of the average position in the synthetic direction (m) after a quench 
on the synthetic tunneling, b. Time dependence of the average lattice momentum (k) along the x direction, c. Schematics 
of the edge-cyclotron orbits, d. Average position in m — x space. The circles in panels a, b, d are experimental data, the 
thin lines connect the points and the thick lines are the theoretical predictions. Experimental parameters: = 27r x 490 Hz, 

t = 27 r X 94 Hz. After the second orbit in panel d, the mismatch between theory and experiment could be possibly ascribed 
to an accumulation of integration error in the data analysis, which amplifies the effects of the assumptions in the model (not 
including, for instance, interactions). 


fall of the edge currents as a function of the Hamiltonian 
parameters [21], and to identify which regimes exhibit 
stronger chiral features. In particular, by varying the tun¬ 
neling rates along x and m, we observe a phase transition 
between a chiral behavior and a non-chiral regime. The 
lattice momentum distribution is measured as a function 
of without affecting other relevant parameters, such 
as Ep and T. Fig. 2d reports the measurement of \J\ 
as a function of Vti/t (circles). As expected, no chirality 
is observed for vanishing when the legs are decou¬ 
pled. Interestingly, chirality is also suppressed for large 
inter-edge coupling Oi ^ t. In the latter regime, the 
largest energy scale in the system is the effective kinetic 
energy along the synthetic direction: this contribution is 
minimized when the fermions occupy the lowest energy 
state on each rung, which does not exhibit any chiral be¬ 
havior. The measured values of \J\ compare well with the 
results of a numerical simulation including thermal fluctu¬ 
ations (shaded area), as described in the Supplementary 
Materials. 

We have then considered a 3-leg ladder, which is the 
minimal configuration where chiral currents at the edges 
can be sharply distinguished from the behavior of the bulk. 
The experimental procedure is analogous to that employed 
for the 2-leg case, adjusting now the Raman parameters 
in such a way as to extend the synthetic coupling to m = 
+3/2, with O 2 — I. 4 IO 1 (see Supplementary Materials). 
Fig. 3a shows experimental data of n{k) and h{k) for each 
of the three legs for Oi = 27r x 620 Hz and t = 27r x 94 


Hz (f^i/t = 6.60). We observe strong chiral currents in 
the upper and lower edge chains, showing values of J 
with opposite sign as in the 2-leg case (J = +0.079(6) 
for m = —5/2 and J = —0.062(4) for m = +3/2). The 
central leg, instead, shows a much reduced asymmetry in 
n{k) (J = 0.018(5)), signaling a suppressed net current 
in the bulk. This is a direct evidence of the existence of 
chiral states propagating along the edges of the system, 
which leave the bulk mostly decoupled from the edges (see 
Fig. 3c). This behavior is akin to what is expected for 
a fermionic system in a Harper-Hofstadter Hamiltonian. 
Bulk states exhibit only local circulations of current, which 
average to zero when all the different states enclosed by the 
Fermi surface are considered. Only the edges of the system 
experience a nonzero current, since there the chiral nature 
of the states prevents that cancellation effect to occur. We 
note that in the ribbon geometry of the experiment the 
bulk reduces to just a single central line. Nevertheless, the 
behavior discussed above is clearly present and detectable 
in the experimental signal. Actually, the small width of 
the ribbon favours the observation of edge states, given 
the large boundary/surface ratio of the system, which is 
reflected in a substantial population of states with edge 
character. 

Fig. 3c shows the values of J for the three different legs 
of the ladder as a function of Vti/t. The results illustrate 
the role of the bulk-edge coupling: similarly to the 2- 
leg case, for small coupling chirality is very weak, and 
increases as Oi/t approaches ^ 3. The theoretical curves 
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show that further increasing Vti/t eventually leads to an 
attenuation of the signal because of the effective coupling 
between the edges, which smoothens the chiral features of 
the system. We observe a substantial agreement between 
experiment and theory for the range of Oi/t that can 
be explored in our setup. The nonzero current in the 
bulk can be ascribed to the different couplings Qi and ^ 2 ^ 
and to a residual light shift which breaks the symmetry 
between the two edges (see Supplementary Materials). 

Finally, we have performed additional quench dynamics 
experiments that directly evidence the properties of chiral 
transport along the edges. A system of lattice fermions is 
prepared in an initial state with zero average momentum 
on the lower m = —5/2 leg of a 3-leg ladder. We then 
perform a quench by suddenly activating the complex tun¬ 
neling in the synthetic direction. Fig. 4a shows the time 
dependence of the average position in the synthetic direc¬ 
tion (m), measured by optical Stern Gerlach detection [20] . 
Fig. 4b shows the time dependence of the average lattice 
momentum {k) along x, measured by time-of-ffight imag¬ 
ing of the whole cloud. Fig. 4d shows an experimental 
reconstruction of the average orbit on the ribbon surface 
as a plot of (m) vs. the average position in real space {x). 
The latter has been determined by evaluating the average 
velocity along x from the knowledge of the energy band 
dispersion vs. lattice momentum, and then performing an 
integration in time (see Supplementary Materials). The 
experimental data are in very good agreement with the 
theoretical predictions, shown by the lines. The dynamics 
displays a strong chiral character, demonstrated by the 
in-phase oscillations in Figs. 4a-b and by the orbits in 
Fig. 4d. Under the effect of the synthetic magnetic field 
the fermions move according to cyclotron-type dynamics, 
which is however naturally truncated by the synthetic 
edge, giving rise to a skipping-type orbit [inilll]. This 
dynamics is effectively damped, already at the theoretical 


level (see Figs. 4a-b), as a result of the average over many 
different fermionic trajectories, which also determines a 
reduction of the average orbit radius to less than one 
lattice site (see Fig. 4d). This is remarkably different 
from the behavior of a non-interacting Bose gas, which 
would occupy a single condensed wavepacket undergoing 
undamped oscillations. 

In conclusion, we have reported the existence of chi¬ 
ral edge states with neutral fermions in a quantum Hall 
ribbon pierced by an artificial gauge field. Our approach 
can be extended to wide ladder systems with up to 2/ + 1 
legs, providing a setting for the investigation of both 
edge and bulk 2D topological matter complementary to 
recent works on Chern insulators [25]. This would allow 
the study of the combined effect of interactions and syn¬ 
thetic gauge fields, a fundamental ingredient for fractional 
quantum Hall physics, in a controlled manner, potentially 
leading to the realization of novel states of matter in 
ladder systems such as, e.g., chiral Mott insulator states. 
Moreover, the fiexibility offered by the present scheme 
allows engineering arbitrary lattice patterns in ladder sys¬ 
tems, including disorder and constriction. This opens the 
door towards the realization of interferometers for chi¬ 
ral liquids, investigating their transport properties, and 
the possibility of implementing interfaces between chiral 
edges potentially hosting exotic non-Abelian anyons such 
as Majorana-like states m- 
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S.I. EXPERIMENTAL SETUP 

The starting point for our experiment is a spin- 

polarized degenerate Fermi gas with Aat — 1.6 x 10^ 
atoms at a temperature T ^ 0.2Tf (where Tp is the 
Fermi temperature). Quantum degeneracy is achieved 
by forced evaporation of a (m = —5/2) -h (m = +5/2) 
nuclear spin mixture in an optical dipole trap with mean 
geometric frequency u/27r 80 Hz. After evaporation, 

the atoms in the m = +5/2 state are removed by a 
resonant laser pulse and a spin-polarized Fermi gas in the 
m = —5/2 state is left. The atoms are then confined in 
a three-dimensional cubic optical lattice with periodicity 
d = Al/ 2 = 380 nm. The three lattice depths are set 
to Vocc = 6.5F^i? and Vby = y^z = 30F^i^ (where = 
h^/2M\\ is the recoil energy, h is the Planck’s constant 
and M is the atomic mass). Along y and z the tunneling 
rates — 1 Hz) are negligible on the timescale of 

the experiment, leading to the realization of an array of 
~ 1000 independent ID fermionic chains characterized 
by a longitudinal harmonic confinement with frequency 
uJxI2t\: ^ 55 Hz. The dynamics is allowed only in the 
shallow lattice along the x direction where t/27r ^ 90 Hz. 
The lattice occupation is less than 1 atom per site on the 
central chain, ensuring that the atoms are in a metallic 
state in the lowest energy band. 


S.II. REALIZATION OF THE FERMIONIC 
HALL RIBBONS 


In this section we describe the laser configuration that 
produces a synthetic magnetic field in the Hall ribbon. 
The complex tunneling along the synthetic dimension 
rh is implemented with two off-resonant A = 556 nm 
laser beams with angular frequencies uj and uj + Aoo. We 
choose a detuning of +1.87 GHz with respect to the 
narrow intercombination transition ^Pi{F' = 7/2) 

in order to reduce the inelastic photon scattering rate. 
The different periodicity of the Raman coupling dR = 
A/[2 sin(6>/2)], with 0 being the relative angle between 
the two beams, with respect to the lattice spacing d, gives 
rise to a non-zero Peierls phase (f = 27r{d/dR). In our 
setup we choose 0 = 19.5° that, considering the projection 
along the longitudinal axis of the fermionic chains, leads 
to a flux ip = 0.37TT (see Fig. SI). The sign of the 
flux can be reversed by swapping the frequencies of the 
two Raman beams. The Raman beams couple up to 
three nuclear spin projections of the ^Sq ground states 
m = (—5/2, —1/2, +3/2) via Am = ±2 processes. These 
states play the role of lattice sites in the extra-dimension 



FIG. SI. Optical setup for the realization of the hybrid 
real+synthetic lattice. Red arrows: optical lattice beams. 
Green arrows: Raman beams. 


m, while the Rabi frequencies of the two-photon processes 
are equivalent to the intersite tunneling amplitudes. The 
atoms are subjected to a Hq = 152 Gauss (real) magnetic 
field along z generating a linear Zeeman splitting Az = 
31.6(7) kHz between adjacent nuclear spin components. 
In order to excite cr+/cr“ Raman transitions (Fig. |S2| ), 
the beams frequency difference is set to Au;/27r 2A^ 

and the polarization is carefully chosen depending on the 
ladder configuration (see below). 

The Raman Rabi frequencies and the light shifts are 
determined by considering the weighted sum of the matrix 
elements over the hyperfine structure F' = (7/2, 5/2, 3/2) 
of the electronic excited state: 
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FIG. S2. Level diagram for the Am = ±2 Raman transitions 
providing the tunneling mechanism in the synthetic dimension 
of the Hall ribbon. 
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where q = (+, —,0) is the photon polarization index, Ap' 
is the single-photon detuning with respect to the F' state 
of the ^Pi manifold and QF is the single-photon Rabi 
frequency related to the i = 1,2 Raman beam and q- 
polarization component. Here a is the synthetic lattice 
site index (see also Fig. 1 of the main text) and 
dicates the coupling between a and a 1. With these 
definitions and considering also the overall harmonic con¬ 
finement along X, the total Hamiltonian of the system 
can be written as: 


N 3 


^=EE[ F h.C.) 


j=l a=l 
N 2 


j=l a=l 
N 3 


E E + Fc.) 


E E + Ar2 E E 


N 3 


j=l a=l 




j=l a=l 
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(S.3) 


In order to take into account the modification of the 
atomic density caused by the interparticle repulsion, we 
consider a slightly lower effective trap strength, Wx — 
0.025N‘^t (corresponding to a trap frequency ^ 40 Hz 
along x). Actually, the results depend only weakly on the 
particular trap frequency used in the calculations. 

In order to calibrate the parameters of the system, we 
set the optical lattice beam along x at Vqx = so that 

the coherent dynamics occurs only along the synthetic 
dimension (t ^ where we observe coherent spin 

oscillations up to 10 periods. Under these conditions the 
kinetic energy and harmonic confinement can be neglected 
and the Hamiltonian can be written as a 3 x 3 matrix in 
a frame rotating at frequency 2A^: 



FIG. S3. Timing of the experimental procedure (see text for 
details). 


elements (^i — h, ^ 2 , ^3 +^) = (0? 0? 2.64)Ui in the rotating 
frame . Under this condition the population of the third 
state is at most a few percent, and can thus be neglected. 

3-leg ladder. In order to resonantly couple all three 
states, we use an out-of-plane linear polarization i = ( 5 + + 
i- +^o)/In this case the light shifts are approximately 
the same for all the states (apart from corrections induced 
by the Zeeman shift of the excited state ^Pi) and by 
setting ^ — ^2 we obtain diagonal matrix elements 

(Cl ~ ^5 C 2 , C 3 + ^) = (O 5 O 5 0.16)Oi. Under this condition 
the three states are all substantially populated. 


S.III. ADIABATIC LOADING PROCEDURE 


H/h = 


[^ 1-5 lli/2 0 \ 

f 2 i /2 ^2 $22/2 , 

0 ^ 2/2 C 3 + 5 ) 


(S.4) 


where h is the detuning with respect to the two-photon 
resonance and U 2 = 1.41Ui. By tuning the polariza¬ 
tion of the two Raman beams we can choose whether to 
implement the 2- or the 3-leg Hall ribbon. 

2-leg ladder. In this case, we choose i = ( 5 + + 
5_)/v^, i.e. linear polarization in the plane xy^ for both 
Raman beams. By setting h = — <^ 2 , we let the first 

two states to be resonantly coupled and exploit the light 
shift to isolate the third one, obtaining diagonal matrix 


In order to adiabatically load the equilibrium state of 
the Hall ribbon we start with a spin-polarized m = — 5/2 
Fermi gas and load it into the optical lattice (Fig. [S^ with 
a 150 ms exponential ramp to the final values V^x = 
and Voy = Vqz = SOEr. During the lattice loading, the 
dipole trap is switched off with a 40 ms exponential ramp 
to decompress the atomic cloud and reduce the overall 
harmonic confinement to uJxI^f 55 Hz. After 5 ms, 
we switch on the Raman beams with an initial detuning 
6in ^ — 25Ui and perform an exponential frequency sweep 
of the form: 


6{t) — 6in + {Sf — Sin) 


1 - 

— p ^adiab I F 


(S.5) 
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FIG. S4. Simulation of the adiabatic loading of the Hall 
ribbon. The solid lines show the calculated lattice momentum 
distribution for the m — —5/2 leg in a th ree-leg system at 
different times while the ramp of Eq. (S.5) is executed. The 
parameters for this evolution are Tadiab — 40 ms and r = 12.5 
ms, while Qi/t = 3.65. The dashed line shows the lattice 
momentum distribution as calculated from the ground state 
of the Hamiltonian with 5 — 5f. 


atomic velocity distribution EaUH], which is measured 
by absorption imaging after 23 ms of ballistic expansion. 

In order to perform a single-site imaging along the 
synthetic direction, we use a sequence of spin-selective 
laser pulses (“blast” pulses), in resonance with different 
components of the narrow intercombination transition 
^5'o ^ = 7/2), to remove atoms in all the spin 

states but one. The sequence is carried out during the 
first 2.5 ms of ballistic expansion. At this time the (real) 
magnetic field is B^ast = 15 Gauss (see Fig. S3), leading 
to a Zeeman shift ^ 50(r/27r) between adjacent spin 
components in the ^Pi manifold (where F = 27r x 180 
kHz is the natural linewidth of the transition). This 
separation allows us to use two opportunely detuned cr+ 
and (7~ beams to remove the unwanted spin population, 
without causing any heating to the selected spin state 
left in the expanding cloud. After ballistic expansion, 
absorption imaging is performed on the dipole-allowed 
^Sq ^Pi (P' = 7/2) transition at 399 nm (with natural 
linewidth F' = 27r x 28.9 MHz) with a resonant pulse 
along z. 


S.V. IMAGE ANALYSIS 


where 6f = — <^ 2 - The ramp duration Tadiab ranges from 

20 to 80 ms depending on the experimental configuration, 
with r ranging from 5 to 25 ms accordingly. The adia- 
baticity of the whole process is verified experimentally by 
reversing the whole procedure to recover a spin-polarized 
Fermi gas. 

In order to assess the validity of our adiabatic load¬ 
ing procedure we also conduct a numerical study. We 
define a time-dependent Hamiltonian where the detuning 
plays the role of a time-dependent chemical potential. 
Then, using the same exponential detuning ramp used 
in the experimental setup, we simulate the evolution of 
the system starting from a spin-polarized Fermi gas in 
a one-dimensional optical lattice. We therefore compare 
the lattice momentum distribution calculated at the end 
of the evolution with the one obtained from t he g round 
state of the Hamiltonian with S = 6f (see Fig. S4). The 
comparison betwe en th e two distributions indicates that 
the ramp in eq. (S.5) creates only a small number of 
excitations in the preparation of the Hall ribbon. 


S.IV. SPIN-SELECTIVE IMAGING 


Lattice momentum distribution and J vs Vt/t, 

The lattice momentum unbalance in a given leg is mea¬ 
sured by taking the average of at least 30 images for each 
spin state with the procedure shown in Fig. |S3| The 
geometric center ko of the lattice momentum distribution 
n{k) is evaluated by acquiring images with Raman beams 
out of resonance, namely with tunneling occurring only 
along the real direction. Indeed, in this configuration 
no chirality is induced in the system, which is purely 
one-dimensional resulting in a momentum distribution 
symmetric around k = 0. All the images are integrated 
along the y direction and shifted by ko in order to obtain 
the n{k) curves for the different legs. In order to remove 
residual gradients or fringes due to imperfections in the 
imaging setup, also background images are acquired, aver¬ 
aged and subtracted from n{k). Finally, we normalize n{k) 
for each spin state in such a way as to have f n{k)dk = 1 
(here, and in the figures of the main text, k is expressed 
in units of the lattice wavenumber kL = 7r/d, so that 
/c = ±1 correspond to the boundaries of the first Brillouin 
zone). The quantity J = h{k)dk for a given spin state 
is calculated by analyzing the average image, while the 
errors are the standard deviation of the statistical distri¬ 
bution obtained after performing bootstrapping on the 
experimental data. 


After the ground state loading, we suddenly switch off 
the Raman coupling, therefore freezing the population 
along the synthetic direction. Then we switch off the real 
lattice with a tmap = 1-35 ms exponential ramp, which is 
slower than the timescale associated to the lattice band- 
gap. This procedure allows us to adiabatically map the 
lattice momentum distribution along the chains onto the 


Edge-cyclotron orbits. In the quench dynamics ex¬ 
periments shown in Fig. 4, in order to study the evolution 
along real space, we determine the average position {x) 
of the atomic cloud assuming the validity of the semiclas- 
sical equation of motion along the real lattice [29]. In 
this framework, considering the lowest band dispersion as 
S{k) = 2t[l — cos{kd)], the velocity of the /^-component 
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of the Fermi sea is 


Vk = T 


1 dS{k) 
h dk 


2td 

h 


sm{kd). 


(S.6) 


Then, knowing the tunneling t, we can measure the av¬ 
erage velocity of the whole cloud at a given time r from 
the experimental lattice momentum distribution n(/c,r): 


{v{t)) = J n{k^r) sm{kd)dk. (S-7) 


Finally, to obtain the average position at a given time r, 
we integrate the average velocity over time: 



{v(T'))dT'. 


(S.8) 


In order to estimate the errors bars of the edge-cyclotron 
orbits in Fig. 4d of the main text, we perform a bootstrap 
analysis over a set of different {x{t)) reconstructed from 
a random sampling of the experimental images. 


S.VI. THEORETICAL DESCRIPTION 


In this section, we provide more details on the simula¬ 
tions presented in the paper. 

J versus Vt/t, In order to quantify the chiral character 
of both 2- and 3-leg ladders, we defined the quantity 

J = [ h{k)dk (S.9) 

Jo 

with 

h{k) =n{k)-n{-k). (S.IO) 


We comp uted this quantity for the Hamiltonian in 
Eq. (S.3) in the central region of the system, with a 
peak density of 0.8/0.9 particles per rung in the center, 
and an approximate size of up to 15 /im (corresponding 
to up to L = 40 lattice rungs with density > 10“^). For 
the data in Figs. 2 and 3 of the main text, we considered 
a single Hall ribbon comprising 75 sites and 20 particles, 
and had the light shifts as the only free parameters in 
the microscopic Hamiltonian for the 3-leg case (in the 
2-leg case the results are not very sensitive to the light 
shift on the —1/2 leg, so that we kept the one extracted 


from the calibration, as described in Section S.H). In the 
3-leg case, we found that the best quantitative agreement 
with the experimental data is obtained by using a set of 
light shifts that slightly differ from the nominal ones. 
In particular, for the simulations displayed in Fig. 3c, 
we used (^ 1 ,^ 2 , Cs) = (0, 0.1, —0.26)fli. The discrepancy 
between these best-fitting values and the nominal ones 
can be easily ascribed to the experimental uncertainty in 
the Raman beams polarization, whose variations have a 
very strong effect on the light shifts. 

In order to take into account thermal fluctuations and 
possible effects of the non-adiabatic loading during the 
lattice ramp, we simulated systems with an average of 
20% thermal excitations above the Fermi sea (both the 2- 
leg and 3-leg case), and performed averages over up to 200 
configurations. These results are approximately recovered 
by investigating a finite-temperature partition function 
with temperature T 0.4t. The width of the theoretical 
curves in Figs. 2d and 3c, representing the effect of 
thermal excitations, has been evaluated as the 68.27% 
percentile interval centered around the mean value of 
the statistical distribution obtained with a bootstrapping 
analysis of 100 realizations of the simulation. 

In addition, we have performed simulations over the 
entire system (up to 2 x 10^ particles in ^ 1000 indepen¬ 
dent ribbons) using a local density approximation in the 
grand canonical ensemble. The results (not displayed) are 
qualitatively the same as for the simple model used here. 

Chiral dynamics at the edge. For the quench dy¬ 
namics illustrated in Fig. 4 of the manuscript, we solved 
numerically the time-evolution of systems at different den¬ 
sities. Indeed, it turned out that assuming a single value 
for the density is usually not sufficient to capture the 
full time-dependent dynamics at the edge. As a matter 
of fact, in the different ladders of the system, each with 
a different density, there is a very different fraction of 
particles participating in the edge dynamics. Whilst not 
affecting drastically the behavior of J, this inhomogeneity 
has strong effects on the single edge-dynamics probed by 
a quench. In order to take into account the inhomogene¬ 
ity effect, we have simulated systems at various system 
sizes up to L = 55 lattice rungs and different densities. 
For the data shown in Fig. 4, we considered a system of 
L = 35 rungs and averaged over a set of different densities 
(0.2, 0.4, 0.6, 0.8). We notice that, when taken singularly, 
none of these density realizations captures the system 
dynamics correctly, while the average gives a very good 
agreement with the experimental results. 
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